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An all-thruster three-axis stabilized attitude control 
system has been designed for the Naval Postgraduate School 
(NPS) satellite bus, ORION. The satellite is a cylinder, 19 
inches in diameter, 35 inches in length, 250 pounds maximum 
mass with 32 pounds for payloads. ORION will be ejected from 
an extended Get-Away-Special (GAS) canister. Launch from any 
GAS can configured expendable booster or the space shuttle is 
assumed. The minimization techniques of Pontryagin have been 
used to derive control laws that support fuel efficent 
operation. A minimum time cost function is applied in the 
acquisition phase to reduce rates to acceptable levels. A 
weighted minimum fuel-time cost function is used during the 
on-station phase.  Bang-Off-Bang control with two switching 
curves is employed outside of a boundary region. Inside the 
boundary region, four-pulse limit cycle control with time 
constants on the order of 100 seconds is applied.  Lifetimes 
in the range of 4 to 28 months are obtained utilizing these 


techniques. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. "While every effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 


additional verification is at the risk of the user. 
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I. INTRODUCTION 


The current ORION design uses spin stabilization, but a 
truly general purpose satellite should provide attitude 
control design options for possible payloads. In surveys of 
likely Navy and Air Force Space Test Program (STP) payloads 
during the past decade, the majority of the requests were for 
three-axis stabilization as shown in Table 1-1 [1]. This 
bears out for the potential customers who have shown an 
interest in the Orion bus. The intent of this thesis is to 
investigate a mass expulsion three-axis stabilization design 
as an ORION bus option. 

A three-axis stabilized Orion bus will undoubtedly be 
nore expensive than a spin stabilized version. Limited 
changes to the original design are desirable. Modified 
deployable solar panels will be necessary to achieve the 
required end-of-life power. Additionally, a revised sensor 
package with more complex sun, horizon sensors and rate gyros 
will be reviewed. This thesis investigates minimum fuel-time 
optimization with respect to thruster operation using 
Pontryagin's Minimization Principle. The effects of external 
disturbances including gravity gradient and aerodynamic drag 
cannot be ignored. Models of a three-axis ORION were 


simulated. 


Table 1-1 SUMMARY OF TYPICAL NAVY/STP PAYLOAD 


REQUIREMENTS [1] 


DataRate 


Orbit 
Inclination 
Instruments 
Control 


32 Ibm 

2.36 ft? 

34 watts 

5000 bits/sec 

200-800 nm circular 
0°-30° or 607-7202 
Particle counter or Lens 
3 AXIS, 2.950 





II. BACKGROUND 


Today most satellites built in Western countries are 
large, complex and multi-mission/payload capable. The cost 
of additional payloads has proven exponential due to many 
factors, of which the most debilitating is that contractors 
are not motivated to organize their program structures 
efficiently. Other factors include: expensive launch from 
only two available sites, the loss incurred when a booster or 
satellite subsystem fails and the increased testing required 
before launch. This puts the free world in a highly 
undesirable position. Reserve assets and the ability to 
produce large numbers quickly i.e., surge capability for these 
high value satellites is extremely limited, if not non- 
existent. In addition, most programs have been stalled at the 
1970-75 technology level due to a general fear of using 
unproven components when the probability of success is already 
questionable and the cost is unacceptable. The concept of 
small, cheap satellites should not work against but with 
current programs. A multi-site launch capability exists from 
ships, submarines, aircraft and mobile platforms. New 
technology can be implemented and space qualified on low value 


vehicles. The ability to build up reserves and to launch 


several satellites at a time becomes viable. The advantages 


are myriad. 


A. CURRENT SMALL SATELLITES 

The present U.S. menu of small proven space platforms 
consists of SPARTAN, NUSAT, GLOMR and LDEF [1]. SPARTAN is 
a retrievable rectangular platform launched from the Space 
Shuttle that can fly alongside for only 2.5 days with a 
payload limit of 1000 pounds. NUSAT (Northern Utah Satellite) 
is also Shuttle deployed but from a Get-Away-Special (GAS) 
canister with a specially designed launch mechanism displayed 


in Figure 2-1. 
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Figure 2-1  Get-Away-Special C 


GLOMR (Global Message Relay) was a DARPA (Defense Advanced 
Research Project Agency) payload which used the NUSAT design. 
LDEF (Long Duration Exposure Facility) is basically a large 
(up to 20,000 pounds), multi-experiment platform that provides 
a facility for extended space exposure. Figures 2-2, 2-3 and 
2-4 illustrate SPARTAN, NUSAT and GLOMR. Though these 
represent advancements in small satellite development, they 
are fundamentally free-flying structures incapable of transfer 


to higher orbits or long duration attitude control [1]. 
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Figure 2-2 Structure for Spartan Spacecraft (1] 
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GLOMR Satellite [1] 


The ability to produce small multi-purpose satellites has 
existed since the dawn of the space age. The need to produce 
a simple yet reliable, low cost satellite bus that is 
available to commercial and military interests alike has never 
been as great as it is now. The Naval Postgraduate School 
(NPS) satellite bus concept not only encompasses propulsion 
and attitude control subsystems, but also standard electric 


power, telemetry, thermal and computer subsystems. 


B. ORION HISTORY 

The history of NPS satellite development began in spring 
1985 when the then Chairman of the NPS Space Systems Academic 
Group, Dr. Allen Fuhs, taught the annual spacecraft design 
course. The project for that particular class of space 
engineering and operations students was the design of a small 
submarine launched satellite [1]. The course was taught two 
more years before Dr. Fuhs retired from NPS. Each of these 
two classes worked on applications of the satellite. One 
student, LCDR Austin Boyd, was so intrigued with the basic 
concept that he pursued a thesis focused on the simple, 
reliable design of three satellite subsystems: structure, 
propulsion and attitude control. LCDR Boyd proposed a Shuttle 
GAS canister deployable satellite, like GLOMR or NUSAT. The 
GAS program is a Shuttle option that allows inexpensive, 
autonomous experiments to fly space available in the cargo 


bay [1]. Due to the unceasing efforts of LCDR Boyd and Dr. 


Fuhs, agencies such as DARPA, SDIO (Strategic Defense 
Initiative Organization) and Navy/Air Force STP (Space Test 
Program) Offices began showing a marked interest in the GAS 
satellite named ORION and in small spacecraft in general. 
This interest prompted the USAF STP office to pursue the 
development of an extended and enlarged GAS canister with an 
improved launch mechanism. The improved GAS can illustrated 
in Figures 2-5 and 2-6 can hold a satellite that is 35 inches 
long, 19 inches in diameter and can weigh up to 250 pounds 
[1]. In addition to the NPS satellite, three other efforts 
were launched in 1986 toward small, general purpose design. 
These included: "T-SAT" by INTRASPACE Company, an advanced 
NUSAT by GLOBESAT and American Rocket Company and an Orion 
type design by TRW [1]. This sudden flurry of activity in the 
small satellite business was capped by an encouraging DARPA 
program initiated in early 1988 called Lightsat. The Lightsat 
program promotes low-cost, lightweight satellites to be 
utilized by operational field commanders [1]. It appears that 
the idea of the small, general purpose satellite has finally 


come of age. 


C. ORION PRELIMINARY DESIGN 

The preliminary design of ORION adopted by LCDR Boyd is 
Simple, reliable and relatively inexpensive with a price tag 
of about $1.5 million. "A general purpose architecture was 


defined as the ability to successfully integrate various 


payloads of the proper size while providing a propulsion, 
attitude control and standardized electrical, data and 


attitude control interfaces." [1] 
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Figure 2-5 USAF/BALL Aerospace Extended GAS Canister [1] 
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Figure 2-6 Comparison of GAS Canister Design [1] 
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General purpose naturally implies that trade-offs must be 
made in the design of the payload in order to reap the 
benefits of a low cost bus.  Trade-offs are weight and size 
driven. The payload of any generic bus is constrained to 
weight and volume limits that do not affect most designs where 
the bus is built around the payload. A brief description of 
ORION subsystems will allow insight into the preliminary 
design. 

The ORION structural subsystem weighs approximately 40 lbm 
and consists of a .62" thick baseplate, three .75" honeycombed 
equipment decks, four extruded 1/16" thick longerons which 
house 78.5" long booms, a 1" thick honeycombed propellant tank 
strongback and four .05" thick semicircular skin panels. The 
baseplate, longerons, booms and skin panels are made of 7075- 
T6 aluminum while the honeycombed items use stainless steel 
facing material. The booms are designed to support 
magnetometers or 2 lbm tip weights to provide stability about 
the longitudinal spin axis. During Shuttle launch, this 
cylindrical structure will experience less than .035" of 
deflection. In compliance with NHB  1700.7A safety 
specifications, the ORION model resonates above 35 Hz at modes 
of 160.8, 178.3 and 244.6 Hz [1]. Figure 2-7 demonstrates the 
basic structural subsystem. 

The propulsion subsystem uses hydrazine monopropellant and 


Shell 405 catalyst. The hydrazine tank is a spherical 
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Figure 2-7 ORION Structural Mock-up [1] 
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positive expulsion type. It is nitrogen pressurized and can 
hold up to 71.5 lbm of monopropellant. |. One 40 Ibf Chrusees 
will be used for orbital maneuvers while six 0.1 lbf thrusters 
will be used for attitude control [1]. Design requirements 
included the ability to transfer from a Shuttle orbit of 135 
nm to 835 nm, the lower limit of the Van Allen radiation 
belts, in addition to a lifetime expectancy of 90 days to 3 
years [1]. 

The ORION power subsystem consists of 1765 in? of ie 
outer skin covered by solar cells and a back-up supply of 24 
NiCad cells rated at 1.25 volts each [1]. An end-of-life DC 
requirement of 60 watts is based on the following assumptions: 
(1) one third of the skin surface area will be exposed to the 
sun, (2) solar cell efficiency of 14%, (3) solar constant of 
125.6 watts/ft^ and, (4) 10% degradation per year. Even with 
these liberal assumptions, end-of-life power occurs after only 
18 months. 

The ORION attitude control subsystem consists of six 0.1 
lbf thrusters, one sun sensor, one horizon sensor and four 
magnetometers mounted on 78.5" booms. Spin stabilization 
about the longitudinal axis with +1 degree of pointing 
accuracy was chosen as the simplest, most cost effective 
attitude control design option for ORION [1]. The fuel budget 
of the attitude control subsystem is limited by transfer 


altitude and effective deployment of the stabilizing booms in 
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this configuration. Figure 2-8 describes the attitude control 


and propulsion subsystems. 
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The remaining ORION subsystems were briefly discussed by 
LCDR Boyd. The thermal control subsystem consists of strip 
heaters and insulating blankets. The health and welfare of 
Orion will be monitored by a computer subsystem and telemetry 
Linke |e Note that the preliminary design assumes Space 
Shuttle deployment which is questionable due to the as yet 
approved use of hydrazine in a GAS can. After all the 
subsystem masses are tallied, a total of about 32 lbm remains 
for payload usage. Figure 2-9 shows the component placement 


for the ORION preliminary design. 
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Figure 2-9 ORION Cross- 
section [1] 
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Table 2-1 ORION COMPONENT MASSES [1] 


Component Mass (lbm) 


Top plate 

Longerons (4) 

Payload Midplate 

Propellant Tank Strongback 
Baseplate 

Launch Restraint Pins(8) 
Pressurant Mid-deck 
Batteries/Boxes (4) 

Earth Sensor 

Propellant Tank 

Hydrazine Control Thrusters (6) 
Orbital Transfer Thruster 
Attitude Control and Payload Computer 
Telemetry Transmitter/Receiver(2) 
Pressurant Gas and Bottle(2) 
Fill and Drain Valve(3) 
Assorted Propulsion Tubing 
Pyrotechnic Valve(2) 

Power Conditioning Circuits(2) 
Payload 

Skin and Fasteners 

Solar Cells 

Booms (4) 

Hydrazine Fuel 
Magnetometers (4) 


Total Mass 
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III. SPECIAL THREE-AXIS CONSIDERATIONS 


A three-axis stabilized ORION will impose additional 
constraints on the structural, propulsion and power subsystem 
design compared to a spin stabilized version. A major 
consideration is the increased fuel usage for a mass expulsion 
attitude control system (ACS). There are several options 
available to extend the ORION lifetime. Although not 
discussed in this thesis, but promoted as avenues of further 
research, there exist other energy efficient three-axis 


attitude control schemes. [2] 


A. ALTERNATE THREE-AXIS SCHEMES 

Momentum bias systems consist of a momentum wheel and 
either mass expulsion thrusters or magnetic torquers to dump 
excess momentum and control motion about the other two axes. 
The momentum wheel supplies gyroscopic stiffness which, in 
essence, cancels unwanted moments about one axis. Significant 
fuel savings result in the thruster moment bias design in 
comparison to the all-thruster system, but the added mass of 
the wheel must be compensated to maintain the same payload 
mass. Magnetic torquers use the torque produced when a 
magnetic dipole cuts through the Earth's magnetic field, but 


they are limited to low Earth orbits. They provide an endless 
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reservoir of control energy as long as current flows in the 
satellite magnets, but they are not preferred if fast reaction 
times are required. The necessity for propellent is 
diminished unless orbital transfer is desired. [2] 

Gravity gradient boom and thrusters or magnetic torquers 
afford a less expensive form of control since the boom is 
passive. Accuracy is reduced compared to the active momentum 
bias system, but fuel is still conserved. 

Reaction wheel systems operate independently about each 
axis just as a three-axis mass expulsion system does. [2] 
Thrusters are employed to dump excess wheel momentum so a fuel 
Savings is realized when equated to the all-thruster system. 
The disadvantages of reaction wheel systems include the 
increased weight, complexity and cost of three wheels which 
could disqualify this system for a small satellite of ORION 


Caliber. 


B. FUEL SAVING METHODS 

In addition to the minimum fuel-time control laws 
discussed and simulated in Chapters V and VI, fuel 
conservation can be practiced by restricting the employment 
of the 40 lbf orbital insertion thruster. Confining usage of 
the 40 lbf thruster to negation of the effects of aerodynamic 
drag could liberate most of the hydrazine for on-station 


@eecicuae control. 
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The ability to transfer between circular orbits is a 
luxury for a small satellite, but AV required is lessened when 
the booster that carries ORION can insert close to the 
requested altitude. The original ORION designers assumed that 
the Space Shuttle would be the only practical launch option. 
In light of the stringent safety requirements for a man-rated 
payload, the shuttle seems to be the least feasible launch 
alternative when dealing with hydrazine monopropellant. Table 
3-1 lists current launch vehicles, payload weights, altitude 


ana) launch costs iG); 


Table 3-1 SMALL SATELLITE BOOSTER OPTIONS 


Booster Name Payload (kg) Altitude(km) Cost(10°) 


Titan 14,030 LEO 
4525 GEO 
Atlas/Centaur 590902275215 LEO(1228.5* 
2330-3145 GTO(i=28.5" 
Delta II 4460 LEO 
1891 GTO 
Scour G=1l 219 555 
*Minuteman III 680 LEO(polar) 
*Conestoga II 724 LEO (1=98°) 
1040 LEO (i=38°) 
*Super Starbird 2080 LEO (10M) 
724 LEO (polar) 
Liberty 1A 7260 LEO 
Pegasus 407 LEO(i=0°) 
271 LEO (polar) 
Foreign Systems 
Ariane IV(France 4200 GTO 
Long March (China) |1267-2489 LEO 
OG: Ne ey: GTO 
*H-2 (Japan) 8978 LEO 
7939 GEO 
Proton (USSR) 20,000 LEO (i=52°) 
2000 GEO 





*not commercially available yet 


20 


If ORION must transfer to a higher orbit, the fuel budget 
will be driven by the AV required. The equation for the 


velocity of a circular orbit is 


V. — (nu/R)* (951) 


where u=GM,=3.986x10” km’/s’, G is the Universal gravitational 
constant and M, is Earth's mass. The radius R=R,+h, R,=6378 km 
is Earth's radius and h is the satellite altitude. A Hohman 
transfer is the most efficient means of transitioning between 


coplanar orbits and is depicted in Figure 3-1. 





Figure 3-1 Hohman Transfer [4] 


The total AV required for a Hohman transfer is 


Dvn AV] + AV, (352) 
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where AV, is the difference between the elliptical transfer 
orbit velocity at perigee and the inner circular orbit 
velocity and AV, is the difference between the transfer orbit 
velocity at apogee and the desired circular orbit velocity. 


The elliptical orbit transfer velocity is 


V = ((24/R) + E)” (33) 


The specific mechanical energy of the transfer orbit is 


E = - p/2a (3.4) 


where a is the length of the semi-major axis of the transfer 


ellipse. The thrust equation for a rocket engine is 


F - mV, * A.(D. - Po) (3259 
m, — mass flow rate 
V, — exhaust velocity 
A, = exhaust cross-section 
p, » exhaust pressure 
p, » ambient pressure 


Specific impulse is a performance measure of the thrust to 


mass flow rate ratio for a rocket 


ISP - F/m (3.6) 


22 


ISP is a function of propellant type. The specific impulse 
of hydrazine is ISP=235 sec. The total impulse is the 


aggregate thrust obtained over time 


I, = f, Fdt = ISP m, (3.7) 


where m, is propellant mass. The total propellant mass 
required to boost a satellite to a higher orbit using a Hohman 
transfer is 


m 2m (1- ee (3.8) 


p O 


where m, is the original satellite mass. The curves in Figure 


3-2 demonstrate the propellant mass required to boost ORION. 


TRANSFER ORBIT FUEL BUDGET 


— 
E: 

— 
Y 
Y 
« 
mk 
a 
d 
E 
% 
= 
o 

e 
e 
E 
A, 


600 1000 1200 1400 1600 


Transfer Orbit (km) 





Figure 3-2 Hydrazine Required by ORION for 
Transfer Orbits 
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ORION will only carry 71.4 lbm of hydrazine so the fuel 


remaining after transfer is alloted for attitude control. [1] 


C. SOLAR PANELS 

Three-axis stabilized satellite power subsystems generally 
include arrays of solar panels that extend after the transfer 
orbit: Panels allow more direct utilization of the sun's 
energy by maximizing the perpendicular rays of received 
sunlight through proper panel rotation. Solar cell output 
voltage is optimized in this configuration. This is the topic 
of another thesis design study and will not be examined in 
detail in this or later chapters. Solar panels will add more 
mass to the ORION structure and thereby add further size 
constraints to other subsystems. Suggested solar array 
designs for a three-axis ORION consist of cylindrical body 
mounted arrays that deploy or accordion style arrays that 


unfold and extend as shown in Figures 3-3 and 3-4. 


To earth 





Figure 3-3 Extended Body 
Mounted Solar Array [5] 
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Figure 3-4 Accordion Solar Arrays [1] 


Solar panels alter the moments of inertia of a satellite 
when deployed. ORION attitude control design assumes that 
I,=1, which greatly simplifies the control problem by causing 
one axis to act independently of the other two in all regimes 
of motion. For I,I,I,, coupling occurs between all axes when 
the angular rate errors are large. An axially symmetric array 
would alleviate the coupling in one axis. 

The additional constraints on fuel and power budgets 
levied by a three-axis satellite are not insurmountable, but 


they must be investigated in more detail in future ORION 


design analyses. 
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IV. THREE-AXIS ATTITUDE DETERMINATION 


Satellite attitude control systems (ACS) consist of 
angular position and rate sensors or estimators and control 
actuators as shown in Figure 4-1. The sensors and estimators 
feed back angle and angle rate information in the form of 
error Signals to be processed as required through a programmed 
control law. The actuators then react to commands from the 
ground or on-board microprocessors to implement corrections 
in attitude. The satellite is affected by unwanted internal 
and external torques such as gravity gradient, aerodynamic 
drag and thruster misalignment. An ACS conforms to the 
equations of motion and reacts to known and unknown 


disturbances. [2] 


Disturbance 
torque 


Attitude angle 


ref - Attitude 
Spacecraft angle 6 
Control law Actuator attitude 
5 dynamics 


command 6 





Figure 4-1 Attitude Control Block Diagram [2] 
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A. SATELLITE DYNAMICS 

In modeling the ORION ACS, rigid-body dynamics will be 
assumed. Rigid body dynamics are used to approximate systems 
that have limited flexibility. The equations of motion for 
a satellite can be described by the movements of a point mass 


in space as shown in Figure 4-2. 





Figure 4-2 Motion of a 
Point Mass [2] 


The mass P can experience rotation or translation with respect 
to the inertial reference frame XYZ. The position vector of 
P is 


R= R,+ 2 (4.1) 


where R, is a vector from the XYZ inertial frame to the xyz 


body frame and r is the position vector of P in the xyz frame. 


r= xi + yj + zk (4,2) 
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The absolute velocity of P is 


ea a (4.3) 


where r includes linear and angular velocity in the body 


frame. 


r=" o er (4.4) 
Substituting (4.4) into (4.3), the absolute velocity becomes 
V -— V. t VL. ToT aE (4.5) 


The absolute acceleration of P is the derivative of absolute 


velocity: 
a= a, t aena + 20 Xx Vi, + ooo ae (4.6) 


In order to simplify future equations, two assumptions 
will be made. The xyz frame is internal to the spacecraft 
body so that V,,-0 and the spacecraft center of mass (CM) 
coincides with the center of the body frame so that f,rdm=0. 
The linear momentum of a rigid-body is defined as the integral 


of the absolute velocity times the differential mass. 
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P - f,(V, * o x r)dm = nV, (4.7) 


Angular momentum of a rigid-body is defined as the moment of 


the linear momentum or the cross product of r and P. 


H, - f,(r x V)dm = fa r x (o x r)dm (4.8) 


The angular velocity vector is equal to 


Qo = w,i t oy,j + w,k (4.9) 


Moments and products of inertia are analogous to the mass 
of a linear system but relate to mass distribution for a 
rotating rigid body. Inertia is the sum of the products of 
the component masses and their distance form the chosen axis 
of rotation. 


ren r. (4.10) 
Moments of inertia become 
= 2 2 
Des OF Gee ae Am (4.11) 


Ben r z) 


I, = = m,(x* + y’) 
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Figure 4-3 shows ORION fuel and dry mass densities along the 


X and Y axes. 
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Figure 4-3 Mass Distribution a)X-axis 
b) Y-axis [1] 
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Table 4-1 lists preliminary subsystem masses and principal 


moments of inertia for ORION. 


Table 4-1 SUBSYSTEM MOMENTS OF INERTIA [12] 


Subsystem 


Structure 
Propulsion 
Telemetry 
Attitude Control 
Computer 


Electrical Power |13.84 


Thermal 
Payload 


Total (kg-m’) 
Total (Slug-ft°) 


Without booms 
Total (kg-m‘) 
Total (Slug-ft^) 


Mass (kg) 


Ix 


2.4652 
1.6218 
0.6089 
9.3621 
0.3482 
07379 
0.0000 
120872 


16:23 
qe 


Iy 


2.4652 
1.6647 
0.5500 
9 3927 
0.3482 
0.5016 
0.0000 
1.1117 


16.02 
11.82 


IZ 


1.4148 
0.8844 
0.0652 


18:6291 


0.0214 
0.3011 
0.0000 
0.0916 


21.04 
15.97 





The products of inertia below are zero when the principal 


moments of inertia are chosen as rotation axes. 


After solving (4.8) 


angular momentum becomes 


and recognizing 


(4.9) 


and 


AS 


(412) 


The inertia matrix is called the inertia tensor. To simplify 
the design of ORION, the inertia tensor will be assumed 
diagonal with zero cross terms. 

The kinetic energy of a rigid body is defined as the 


integral of the dot product of the absolute velocity. 
T - if, V.V dm (4.13) 


This thesis will focus on rotation only. Substituting (4.5) 


into (4.13) the kinetic energy reduces to 


o, (4.14) 
= 


(4.15) 


A moment about the satellite CM is the cross product of 


distance and force. 


Fa 
I 


- = H + f,R x rdm- H, (4.16) 


= Herel tox H. 
Angular momentum in the body frame can be expressed as 
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Me 12 + 4) eek (4.07) 


where xyz coincides with the principal moments of inertia. 


From (4.12), the elements of the angular momentum vector are 


H, = Lo, (4.18) 
H, = Iyoy 
HI 


After some manipulation, (4.16) is transformed into the Euler 


Moment equations. 


FA 
x 
| 


A o E (4.19) 


Iw 


y2y + w,0, (I, - I) 


J 
N 


M, = I w, + Gray (Ly =. 19) 


Equation (4.19) will be used extensively through out this 
thesis to simulate ORION attitude control. 

An inertial coordinate system with X,Y,2, fixed at the 
center of the Earth as shown in Figure 4-4 is used as a 
reference frame when describing satellite motion. The roll 
axis, X, describes the orbital velocity vector with magnitude, 
wo, which is assumed constant for a nearly circular orbit. 


The yaw axis, Z, extends from the satellite CM toward the 
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Figure 4-4 Spacecraft Coordinate System [2] 


Earth CM. The pitch axis, Y, is perpendicular to the orbital 
or XY plane. Errors caused by disturbance torques are 
described by the angles 4, y and @ about the roll, yaw and 
pitch axes, respectively. In unperturbed form, the Euler 


moment equations become 


lo; D ME (4.20) 


Eo; (I. EI) 


lzo, (I I,) yo, 
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The right hand side of (4.20) exhibits the coupling that 
occurs when angular error rates are large. For on orbit 
control of a three-axis stabilized system the angular 


velocities are assumed small which reduces (4.20) to 


Iv, = O (4.21) 
Io — 18 
Iv, = O 


Satellite acquisition is a stage in which initial large 
attitude adjustments are made to correct for booster 
separation. In this phase the satellite frequently 
experiences a tumbling motion which translates to large errors 
in angular position and rate. Excess rates must be negated 
before the on orbit control law can be implemented. Angular 
velocity errors as viewed from the inertial frame can be 


expressed as 


(4.22) 


€ 
E 
ll 
De 
I 
€ 
O 
= 


€ 
N 
ll 
e 
+ 
€ 
o 
D- 


Replacing rate errors in the Euler moment equations yields 
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= 
x 
| 


- L($ - o9) - (Ll, DA (4.23) 
M, = I, 


M, = IL(v 9 0,9) + (I, I)o E RD 


B. EXTERNAL DISTURBANCES 

The major external disturbances for low altitude orbits 
are caused by the Earth's gravitational field and aerodynamic 
drag. Solar pressure contributions are considerd negligible 
at low altitude but a dominant factor beyond geosynchronous 
altitudes. 

1. Gravity Gradient Torque 

Gravity gradient torques can be used as a passive form 

of three-axis stabilization for Earth pointing payloads if a 
boom is attached to the satellite. This torque results from 
the interaction between the earth-satellite masses. [2] The 


gravitational moment can be expressed as 


M = fa r x F dm (4.24) 


where the gravitational force is 





B (OR, T = a e (4.25) 
ps 


Substituting (4.24) into (4.23) results in 
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3 
M; = —, f(r x Re)(r - R,) 
R 


e 


(4.26) 


where R,, the vector between Earth-satellite centers of mass 
as shown in Figure 4-5 is expressed in the satellite body 


frame as 


R. = R.(-sinfi + singcosg) + cos¢cosék) (4.27) 





Figure 4-5 Gravity Gradient Moment on a 
Satellite [2] 
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For small angle corrections equation (4.26) reduces to 


Lv 
Ms ~ EC) (4.28) 
0 


2. Aerodynamic Drag Torque 

Aerodynamic drag is caused by the friction between the 
satellite as it orbits at an angular rate of w, and the 
Earth's atmosphere. It is felt most below 500 km. The 
effects of drag cause a general decrease in angular velocity 
thereby resulting in decreasing altitude or orbital decay. 
The aerodynamic drag force is a function of air density, p; 
velocity, V; surface area perpendicular to the flight path, 
A; and a dimensionless number called the coefficient of drag, 


Cp. [8] 
Ep = Ion ER (4.29) 


The orbital lifetime of a satellite is a function of altitudes 


h, and eccentricity, e. [7] 


3/2 
60378 Fh (4.30) 
je 


L = NP = 1.15185x10" N | 
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where N is the estimated number of revolutions and P is the 


The values of N can be obtained directly 


period of the orbit. 


from Figure 4-6 where the ballistic coefficient for ORION is 


Note that the eccentricity for a 


m/(CjA)2132.87 kg/m. 


orbital lifetime is 


circular orbit is zero, 


For ORION, 


e=0. 


The linear 


shown in Figure 4-7 as a function of altitude. 


velocity of a satellite equals 


(4.31) 


wR = (u/R)* 


ll 
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ENS NAS Ls 
MITIANOSNAN NONE I N 
MIT IAN NUS, III] 
A RSS HTA" 
[SO 
ES TM 
CRUISER LT HTL, 


“pa. e i = x N ¿oo » N O70 FT N COMO T N 
”=O00 © eo 


((¿0/54)1/SNOILNTOA3H) 
1N312133302 2I1SI11v8/3WI 13311 31v WIXOHddVv 


z 
= 
= 
Ps 
2 
ev 
L 
Pr. 
ev 
a 
c 
uJ 
a 





Orbital Decay [7] 


Figure 4-6 
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Figure 4-7 ORION Lifetime due to 
Aerodynamic Drag 


where R = R, + h, the altitude above Earth CM. The drag 
coefficient is between 1 and 2 in SI units. The value used 
for ORION will be C2. The surface area of ORION 


perpendicular to the flight path is A=.429m*. This reduces 
(4.29) to 


Fy = 1.71 x 107R (4.32) 


ORION drag torque will can be summed up as 


Taras DASS (4.33) 
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where L=.076m is the ORION center of pressure location along 
the y axis.[8] The effects of aerodynamic drag and 
gravitational gradient are simulated for a three-axis ORION 


in Chapter V. 


C. SENSORS AND ACTUATORS 

The ACS sensor suite is dictated by the type of 
stabilization, pointing accuracy, altitude, orbital parameters 
and cost constraints of the satellite budget. There are 
potentially three distinct stages of attitude control that 
ORION must handle. These include transfer orbit, acquisition 
and small error correction which includes altitude 
maintenance. If a transfer orbit is required, the satellite 
Will most likely be spinning during this phase so different 
types of attitude determination sensors will be needed to 
accommodate spinning and  non-spinning modes. In the 
acquisition stage, the satellite experiences large error rates 
which must be negated before accurate pointing can be 
sustained. The small error stage is the steady state 
operation that the satellite will undergo for the majority of 
its active life. 

The attitude of a spacecraft in inertial space can be 
completely realized by two vectors as long as they are not 
colinear. Point sources such as the sun and stars afford 
excellent reference vectors while the Earth and planets 


present extended sources. Pitch, roll and yaw attitude errors 
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relative to the desired satellite orientation must be 
determined in addition to error rates about these three axes. 
Combinations of horizon and sun sensors can furnish position 
error while rate sensing gyros or error estimation routines 
can be used to resolve rate errors. 

The current selection of attitude sensors for a spinning 
ORION consists of sun Sensors, horizon sensors and 
magnetometers. The sun and horizon sensors provide the 
primary angular error output while the magnetometers act as 
back-up sensors. The spinning motion of the original 
satellite simplifies the number and cost of primary sensors 
by producing a periodic search motion. In a three-axis 
stabilized system this scanning motion must be integrated into 
the instrument in the case of the horizon sensor while 
multiple sun sensors with wider fields of view must be 
employed. 

1. Horizon Sensors 
Earth or horizon sensors are infrared devices that 
detect the 14-16 um radiance profile of the earth as the space 
to earth and earth to space transitions occur. Horizon 
sensors can be used to ascertain pitch and roll errors. [2] 
Scanning or radiation balance detectors are employed in three- 
axis designs to provide a sweeping action across the disk of 


the Earth. Scanning sensors use two beams centered on a 
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reference line to scan the Northern and Southern hemispheres 


in an east-west direction. The scan pulse width is 


pw- d/(2hw,) (4.34) 


where d is the distance between horizon crossings, h is the 
satellite altitude and w, is the scanning rate of the sensor. 
The difference between in-crossing to reference and reference 
to out-crossing pulse widths for each beam can be used to 
determine pitch error while the difference between Northern 
and Southern scans is proportional to roll error. An 
incremental pulse generator outputs a pulse at fractions of 
a degree for the known scan rate. Pulse width is computed by 
adding pulses from acquisition of signal (AOS) to reference 
then from reference to loss of signal (LOS). [7] Typical 


output from a scanning horizon sensor is depicted in Figure 
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Figure 4-8 Output from a 
Scanning Horizon Sensor [7] 
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2. Sun Sensors 

Since two non-coincidental directions are needed to 
measure alignment uniquely, it is necessary to find another 
reference vector besides that from the Earth. The sun or 
another star are logical alternatives since they can be 
approximated as point sources. Stars are the best inertial 
reference points and are not coplanar unless the inclination 
is zero and the orbital path is in the equatorial plane. Star 
sensors are orders of magnitude more complex, expensive and 
bulky than sun sensors which owe their simplicity to the close 
proximity of the most radiant body in the solar system. 

Sun sensors are the most commonly used attitude sensor 
due to their relatively low cost, lightweight and broad 
application. [2] They can render accurate position error 
updates in harmony with horizon sensors when the satellite to 
Earth and satellite to sun vectors are non-parallel. Vector 
alignment can occur up to twice an orbit at local noon and 
midnight. Sun sensors are basically silicon solar cells with 
specialized masks for varying functions. Analog and digital 
sensors are available to measure the angle between the sun 
line and the pitch axis or sun presence in a particular field 
of view for instrument protection. The sun line vector is 


defined as 


S = sinecosól +sin$J cosacoséK (4.39) 
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where $ is the sun declination and a is the orbital angle of 
the satellite measured form local noon. Yaw errors about the 


pitch axis in the inertial frame can be expressed using 


j = -pI + J + 4K (4.36) 


for small angles. The sun angle, 9,, is obtained by taking 


the dot product of the sun vector and j 


cosé, = S - j = -y¥sinacoss + siné + ¿cosacos$ (d) 


As the yaw error changes the sun angle varies 


acos? ., 


àv 


= -sinacoss 
j (4.38) 


thus yaw error measurement is obtained. [2] 
3. Magnetometers 

Magnetometers are not a required sensor for current 
or future ORION designs. The original design utilizes booms 
with 2 lbm weights at the end to enhance the stability of the 
spinning satellite by altering the principal moments of 
inertia. Three of the tip weights are magnetometers which 
will provide a back-up attitude determination capability that 


dead weights cannot. Magnetometers are solenoid coils that 
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obey Faraday's Laws of Magnetic Induction. When placed in a 
time-varying magnetic field, solenoid coils will produce a 


voltage 


V=$ E. dl = -AN¿(aB/dt) (4.39) 


where E is the electromotive force, A is the cross-sectional 
area of the solenoid core, N is the magnetic permeability of 
the core and B, is the magnetic field component along the 
solenoid axis. The voltage is proportional to the magnetic 
field magnitude and direction. On-board control using 
magnetometers requires knowledge of the Earth's magnetic field 
at the orbit altitude. Knowledge of the magnetic field is not 
complete at low altitudes and too weak above 1000 km since the 
field strength decreases as the cube of the distance from the 
barch*s CM. [7] 
4. Rate Gyroscopes 

In the acquisition phase, some means of determining 
error rates is required before they can be cancelled. Rate 
sensing gyroscopes, Luenberger observers or Kalman filters for 
the sun and horizon sensor output are viable solutions. Rate 
integrating gyros (RIG) in triplets are the most accurate 
instruments used to detect error rates and positions about 
three axes. The output of a RIG is the current required to 


maintain the gyro gimbals at a null or zero position. The 
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current is proportional to the torque needed to zero the 


rates. The total angular momentum of a gyro is 


H = LS + 1,90 (4.40) 


where L is the magnitude of the angular momentum of the rotor 
about the spin axis and I,0 is the magnitude of the angular 
momentum of the gimbal system about the output axis. The 


nulling torque is 

T, = (dH/dt),,, * e x H (ae) 
Angular velocity in the gyro inertial frame is 

w = wI + w0 + wS (O) 


where I is the gyro input axis, S is the rotor spin axis and 
O is the output axis. Figure 4-9 depicts gyroscope operation 
and the physical location of the axes described by equations 
(4.40) and (4.42). The differential equation of motion for 


a gyroscope is 
1,9 + DO + K9 - vıL = 0 (4.43) 
where D is a viscous damping term and K is a restoring term. 


[7] 
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Figure 4-9 Single Degree of Freedom 
Gyroscope [7] 





5. Hydrazine Thrusters 

The functions of station keeping and attitude error 
correction are performed solely by control actuators called 
thrusters or reaction jets in this design. The thrusters on 
ORION use the monopropellant hydrazine with the Shell 405 
catalyst for these functions. The hydrazine storage tank is 
pressurized with nitrogen. A basic hydrazine reaction jet is 
pictured in Figure 4-10. 

The lifetime of a three-axis stabilized satellite is 
more acutely limited by the operational usefulness of the 
Shell 405 catalyst than in a spin stabilized vehicle where 
fuel usage can be more accurately predicted. The reaction 
jets on a three-axis stabilized satellite will be subjected 


to intermittant firing or low duty cycle firing.  Thruster 
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Figure 4-10 Basic Thruster Design [2] 


performance is largely determined by the temperature of the 
catalyst bed. During steady state operation, intermittant 
firing of the jets will not allow the catalyst bed to maintain 
sufficient warmth to facilitate the most efficient thrust 
profile between firings. Typical values of specific impulse 
for hydrazine during transfer and steady state operation are 
235 and 135 seconds, respectively. The specific impulse is 
seriously degraded in the steady state regime. Ultimately, 
this leads to the conclusion that the catalyst bed in a low 
duty cycle environment must be heated which requires 
considerably more power than just that used for opening and 
closing valves and heating the hydrazine to prevent freezing. 
[2] 

In addition to the catalyst bed temperature, thruster 
performance is affected by impurities in the hydrazine that 
accumulate in the catayst bed after each firing and variation 


in the nitrogen pressure as the propellant is depleted. [2] 
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Figure 4-11 shows the linear relationship between thrust 


obtainable and pressure. 
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Figure 4-11 Thrust versus Tank 
Pressure [2] 


A typical thruster profile is illustrated in figure 
4-12. This profile assumes that physical conditions are 
acceptable, but that there is a finite delay between the time 
that the thruster is commanded on and the time that it reaches 
maximum thrust and likewise for the command off sequence. The 
time from t, to t, and t, to t, is a function of the valve 
opening and closing action. In Chapter V, the thruster 
profile is modeled as a pulse function for simplicity. 

The current thruster selection for ORION includes a 
40 lbf jet for orbital insertion and six .1 lbf attitude 


control jets. Figure 4-13 lists design specifications for 


these thrusters. 
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Figure 4-12 Typical Thrust Profile [8] 
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Design Characteristics 


Propellant Hydrazine | 
Catalyst Shell 405/LCH 202 
[ 
Í 
a 
| 


Thrust, Steady State (Ibf) 40— 11 
Feed Pressure (psia) 450 — 100 
Chamber Pressure (psia) 223 — 61 
Expansion Ratio 203 
Flow Rate (Ibm/sec) 0 17 — 0.054 
Valve Wright Components Single Seat 
Valve Power 50 Walts @ 28 vde and 75°F 
Weight (Ibm) 1.95 

Engine 1.35 

Valve 0.60 


Ceara a 


Demonstrated Performance 


Specific Impulse (Ibf-sec/Ibrn) 235 — 203 
Total Impulse (Ibf-sec) 4,350 
Total Pulses 2.298 
Minimum Impulse Bit (Ibf-sec) 0.71 @ 450 psia & 20 ms ON 
Steady-State Firing (sec) 97 








Figure 4-13a Rocket Research 40 lbf Thruster - 
Model 107B [1] 
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Design Characteristics 
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Feed Pressure (psia) 
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Flow Rate (lbmi/sec) 
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Figure 4-13b Rocket Resear 


Model 103C [1] 
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Hydrazine 

Shell 405 

0.252—0.042 

.420—70 

.370—60 

. 100: 1 

0.001—-0.0002 

Wright Components Dual Seat 


9 Watts Max. at 28 vdc and 45°F 


Ores 
220-25 
0.45 


227—206 

33,025 

e 410,000 

0.001 E 150 psia 4 8 ms ON 
64,800 


.1 lbf Thruster 





V. MINIMUM FUEL-TIME PROBLEM 


Three-axis stabilization is inherently more expensive than 
spin stabilization due to the increase in complexity of 
attitude sensors and the energy necessary to maintain the 
desired pointing accuracy. For Orion this translates to more 
fuel needed once on-station than the original spinning 
version. This limits the total AV available for orbit 
transfer and the active on-station control life of the 
satellite. Therefore, the effective operational lifetime of 
the payload could be dictated by fuel constraints before other 
considerations. In order to conserve fuel and provide maximum 
operational lifespan, a minimum control effort solution using 
Pontryagin's minimum principle is explored in this section. 


[3] 


A. ON-STATION CONTROL LAW 

For small angular rate and pointing errors, examination 
of one axis applies to all three axes. Thrusters are the sole 
control device for this design study and act independently in 
this case. The simplified differential equation for one axis 


is 


(5.1) 
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where T, is the thruster control torque and I is the inertia 
tensor with the prođucts of inertia equal to zero. The 


normalized state equations for one axis are 


le) rl aa 0 (5.2) 
= D mac 
X(t) 0 0 x,(t) 1 

where x, is angular position, x, is angular velocity and u is 


the control torque provided by the thrusters. The minimum fuel 


performance index or cost function is 
J(u) = f, (3 *A[u|)dt (5.3) 


where A is a weighting factor applied to the control. As 1-0, 
the performance index becomes the minimum time solution and 
as A2o, the performance index approaches the pure minimum fuel 
solution. Since any minimum effort problem is a trade off 
between elapsed time and consumed fuel both should be 
considered and are represented by (5.3). The performance 
measure must be minimized to obtain the optimal solution for 
given initial conditions and final states, x(t;)-0. Every 
value of the weighting factor provides a minimal trajectory. 
The selection of the weighting factor is determined by 
reasonable transient time and fuel expenditure. The 


normalized admissible controls for thruster operation are 
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juty TITS (5.4) 


where T, is the maximum torque available from the reaction 
jets. The cost function for thruster operation is called 
'Bang-Off-Bang' since the satellite drifts before a final 
control torque is applied to bring it to the desired final 


attitude. The Hamiltonian for (5.3) is defined as 


H(x,u,p,t) = 1+A|ul/+p,x.+p,u (5.5) 


where p(t) are Lagrange multipliers. The costate equations are 


p,(t) ="="0H/ ox, = 70 (5.6) 
p(t) = - 3H/3Xx, = -p,(t) 

which have solutions 
p (t) = C; (5.7) 
p(t) = -cıt + C; 


The Hamiltonian must be minimized on an optimal trajectory, 


so the form of optimal control is 
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u(t) = 0 IS Dees A 
-N P2 > A 


where switching occurs at |p;|»4. The effect that any known 
or unknown torque has on the satellite can be best illustrated 
in the phase plane. This is a plot of angular position versus 
velocity as a function of applied torque. There are 
essentially two switching curves for the given minimum 
trajectory. These correspond to a curve that is dependent on 
the weighting factor, A, and the minimum time curve as shown 


in the phase plane representation in Figure 5-1. 
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Figure 5-1 Double Switching Curve Control 
Law 


If the initial conditions place the satellite in region 1, a 


negative torque is applied until t, when region 3 is entered 


5 


and control is turned off. Positive control is applied at e 
when the satelite reaches region 2, which is the switching 
curve for the minimum time solution. The reverse is true if 
the initial conditions place the satellite in region 2. At 


t, the control goes off, u(t,)=0, and (5.6) becomes 

P2(t1) = A = -cıt, + C? (5.9) 
The position continues to change as a function of velocity at 
turn-Oorft. After t,, control is turned on, u(t,)=0, and the 
satellite ideally moves to the desired position. At t, (5.9) 
becomes 

p2 (t2) = -å = -cıt + C? (5.10) 


Subtracting (5.10) from (5.9) 


(t, - tj) = 2A/c, (5.11) 


Assuming that velocity is constant at turn-off, x,=0 
X;(ti) = X2(t2) = X2(t) 


Solving the Hamiltonian at the switching times and 


substituting (5.6) for p, yields 
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H(ti) 


H(t,) = 0 





H(t,) = 1 - Cc,x2(t,)=0 
= i (5.12) 
X2 (t1) 


wpet itüting (5.11) into (5.12) 


(t? - ti) = 21X) (5.13) 


Solving the state equation, x,(t)=x,(t) 
X1(t;) - X,(t,) = X2(t,) [t,-t,] (5.14) 


From t, to t; the satellite follows a minimum time trajectory 


where 


2 
x(t) = -} A (5.15) 


ExEStrtuting (5.13) and (5.15) into (5.14) 


M = 0 E 


2 [N| 


which is the control law for this minimum fuel problem. 
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Equation (5.16) becomes usable in all quadrants when 
x,=x,|x,| which takes the sign of x, into account. Figure 5- 
2 demonstrates minimum trajectories for various values of a. 
The initial conditions x,(0)=1,x,(0)=0 were used to test for 
minimization of the selected cost function. As shown in 
figure 5-3, the minimum value of the chosen cost function 
occurs at \=2.25 when this value is substituted into equation 


5.16). ! 
) MINIMAL TRAJECTORIES 
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Figure 5-2 Parabolic Switching Curves for 
1=0,1,10,100 


Figure 5-4a is a time response of angular position and rate. 
Figure 5-4b is the time response of control torque. Notice 
that thruster chatter occurs at the end of the cycle, since 


fine control near the desired attitude is not implemented. 
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Figure 5-3 Minimum Fuel-Time Cost Function 
IMA + 2.25|u| jdt 
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Figure 5-4 Time Response of States and 
Control 
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Figure 5-4b shows the "bang-deadzone-bang" behavior of the 
control law. Figure 5-5 is a phase plane representation of 


the same cycle. 
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Figure 5-5 Phase Plane for \=2.25 


B. LIMIT CYCLE CONTROL 

Ideally, the last thrust period will take the satellite 
to the origin, but this is a function of the integration step 
size. A similiar analogy exists for the thruster in the form 
of the minimum impulse bit available. Physically angular 
rates are never zeroed, but they can be maintained in a 
region. Satellite motion is then dictated by a finite limit 
cycle even in the absence of external disturbances. 

While in the limit cycle attitude errors are reduced below 


specified values. The fine tuning region is determined by 
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pointing accuracy requirements, but it is also a function of 
minimum jet impulse for angular velocity and the precision or 
noise contribution of the attitude sensor instrumentation. 
Minimum jet impulse is a measure of the smallest thrust unit 
available when the control valves are immediately opened and 
closed, i.e., the shortest single jet pulse. In order to 
increase the coasting period in the limit deadzone, the 
satellite rate must be reduced to the lowest value possible 
without causing thruster chatter. Chatter is a constant 
cycling of opposing thrusters when limits are exceeded after 
each firing. The limit cycle for the cost function simulated 


in this example is depicted in Figure 5-6. 
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Figure 5-6 Limit Cycle 
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A third switching curve was implemented in the boundary 
region only to minimize fuel expended while cycling. The 
third switching curve has a flatter slope, A=250, and thus a 
wider deadzone. The best deadzone obtainable is again limited 
by the minimum jet impulse available from the thrusters. 
Figure 5-7 utilizes the third switching curve and can be 
compared to Figure 5-6 where only two switching curves were 


realized to demonstrate the improved fuel consumption. 
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Figure 5-7 Upgraded Limit Cycle 


Figure 5-8 shows the time responses of thruster operation 
during fine tuning in both cases. Note the single impulse 
firing that occurs at the edge of the limit boundaries when 


the third switching curve is implemented. The appendix 
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provides a description of the Tutsim program that simulates 


the independent control law about a single axis. 
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Figure 5-8 Fine Control 


C. ACQUISITION PHASE CONTROL 

The acquisition phase is the most critical aspect of 
satellite operation after a successful launch. During 
acquisition, the satellite may be tumbling. Tumbling 
translates to high error rates which will cause coupling about 
at least two axes. In this phase, a minimum time solution is 
desired to cancel the excess rates quickly. [6] The Euler 


moment equations of motion are 
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I,w, = (I5 wes) Oyo, a Ta (5.22 
Lo, = "E, = ee eee 
lo, — UI. AIME E LM 


where T,, T, and T, are the control torques; the produc 
thrust and moment arm. If a single axis of symmetry is 
assumed, as is the case for a cylinder, two of the moments of 


inertia will be equal 


The moment equations become 


€ eœ 
| 


<= ao O E (5.183 
en 
Q, = -aw,w, t U, 


Where a = (I-I,)/I and u;=T;/I. 
Motion about the axis of symmetry is independent of the 
other two axes and is controlled at all times as discussed in 


the previous section. Assuming that o, is constant when u, is 


y 


zero during rate reduction, the state equations become 


“S . j bee : l (5.19) 
: = u(t) 
x, (t) O x,(t) o i 
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where B=aw,, X,=0,, X¿=0,, U¡=4u, and u,=u,. The minimum time 
cost function is 


dE (522.0) 


The Hamiltonian for the coupled system using both controllers, 


u, and u,, is 
H(x,u,p,t) = 1 + PiX - P2X; + pu, + pu, (95721) 


for B=1. The costate equations are 


pı(t) nPE) (5.22) 
PAN ee Hp pi) 
which have solutions 
Pi(t) = cz; cos(tté) (5:723) 
pt Nc Eum ED 


To minimize the Hamiltonian on an optimal trajectory, the form 


of optimal control is 


-sign(p,) 522) 


u, (t) 


u,(t) -sign(p;) 
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The control torque of each set of thrusters will reverse every 
half cycle while the control mix will change every quarter 


cycle. The minimum time switching curves can be derived using 


ax, _X + u, (5.25) 
ax, — Xi zt U, 





$ (-x,+uz) dx, = $ (x2+u,) dx, 
2 2 
ZI N m 0 
2 


Adding u,'+u,? to both sides of (5.25) and multiplying by two, 
yields 
2 2 


(x; u0,)^ £006 Fu) AE (5.26) 


For normalized control torque, u,=+1, u,=t1, (5.26) becomes 
(x, fF I Fa D (5.27) 
Equation (5.26) describes four circles with their centers at 
(+1,+1) in each quadrant. The switching curves are defined 
as the quarter of each circle that intersects the abscissa and 


ordinate in each quadrant as shown in Figure 5-9. Beyond 


these scalloped sections the switching curves are mirror 
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images of the half sinusoidal sections along each axis, but 
minimal error is experienced if switching is implemented about 


the axes themselves. 


A — 


DOUBLE CONTROL SWITCHING CURVES 





Figure 5-9 Double Controller 
Minimum Time Switching Curves 


A phase plane representation is used to illustrate rate 
reduction where w, is the abscissa and w, is the ordinate. 
Figure 5-10 depicts the implementation of the four quarter 
section switching curves in a circular region with a radius 
of two for initial rate conditions of (#1.5,+1.5). Figure 5- 
11 portrays switching about the coordinate axes in the region 
beyond a radius two. In the outer region, u,=-sign(w,) and 
u,=-Sign(w,). Figure 5-12 shows the control torque time 


Response for u, and u,. 
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Figure 5-10 Inner Region Switching 
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Figure 5-11 Outer Region Switching 
a)w,=w,=6 b)w,=-7 w,=5 
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Figure 5-12 Outer Region Control 


Since the system is coupled, only one controller is needed 


to reduce the rates about both axes. If u,=0, equation 


becomes 


X] 


N 


2 
= UX = 222 


N 
N 


Adding u, to both sides and multiplying by two 


which becomes 


2 2 2 
(x; - uj) * x; = Ur 


(x EE A 
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(5.25) 


(5.28) 


(5.29) 


(5.30) 


for u,=+1. Conversely, for u,=0 equation (5.30) becomes 


x, + (x, E 3) (53) 


for u,=+1. Equations (5.30) and (5.31) describe two symmetric 
circles with centers at x,=+1 and x,=+1, respectively. The 
switching curves for the minimum time solution are the half 
circles shown in Figure 5-13. Figure 5-14 demonstrates rate 
reduction for a single controller system. For initial rates 
in the outer region, switching about the axes was implemented. 
Control torque values of u,=sign(w,) for (5.30) and u,=sign(w,) 
for (5.31) were applied. Figure 5.15 demonstrates large rate 


error reduction for both single controller cases: 
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Figure 5-13 Single 
Controller Switching 
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Figure 5-14 Single 
Control for 
w,=0,=+1.5 
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Figure 5-15 Outer 
Region w,=w,=6 
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Figure 5-16 illustrates the control torque time responses in 
the outer region. Comparison of Figures 5-12 with 5-16 shows 
that the double controller system takes half the time using 
the same amount of control effort as the single controller 


system. 
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Figure 5-16 Outer Region Control Torque 


D. LIMIT CYCLE CONTROL FOR ORION 

ORION specifications can easily be implemented in the 
simulation described in previous sections. A three-axis 
stabilized ORION will be oriented such that the yaw axis is 
colinear with the longitudinal axis. The symmetry of the 
principal moments of inertia becomes I,=I, vice I,=I, as in the 
spin stabilized configuration. ORION principal moments of 


inertia from Table 4-1 are 
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2 
„= I, = 6.753 kg-m 


H 
t 


2.776 kg-m 


H 
N 
ll 


Attitude control thruster data is obtained from Figure 4- 
14b. The thrust available from the attitude control reaction 
jets is between .24 and 1.068 Newtons. The minimum impulse 
bit is .0223 N-sec for .025 seconds. The total pulses 
available as a function of the catalyst effective lifetime are 
161,000. The mass flow rate of the thrusters is between .104 
and .453 grams/sec. 

Torque is the product of thrust and moment arm. The 
length and radius of the ORION cylinder are .889 and .2414 
meters, respectively.  Thrusters operate in pairs to prevent 
translational movement so the average torque produced for 


ORION is 


T, = Ty = 2(.445)(.889) = .7912 N-m 


le 
I 


2(.445)(.2414) = .2148 N-m 


The control accelerations are 


A eo 


- 
I 


0774 s 


E 
ll 
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Pointing accuracies of +1° and +.05°/sec were used in the 
example illustrated to define the limit region. The minimum 
fuel-time switching curves were implemented with the values 
of A previously used, A=2.25 and 250 for the outer region 
deadzone and limit region deadzone switching curves, 
respectively. Figures 5-17 and 5-18 illustrate the minimum 
fuel-time control law for yaw and roll-pitch when ORION values 


for torque and moments of inertia are used. 
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Figure 5-17 Yaw Limit Cycle for Limits of 
+] ° IS AO 


Figure 5-19 shows the yaw thruster operation for initial 


conditions of angle=.25 rad and zero rate. 
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Figure 5-18 Roll-Pitch Limit Cycle for 
Limits of +1”, +.05°/sec 
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Figure 5-19 Yaw Thruster Operation 
angle=.25,zero rate 
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It is difficult to predict fuel usage up to the 
operational phase of flight for three-axis control. The 
approximate limit cycle lifetime with all thrusters operating 
is 

Total Fuel Available 


I, — _____ > At (5.399 
Fuel per Control Cycle 


where At - one limit cycle. [1] In the limit cycle, with no 
environmental disturbances added, individual pairs of 
thrusters pulse on four times per cycle. Using a mass flow 
rate of .104 grams/s for .025 sec, the minimum pulse time, 
gives a total fuel per control cycle of between .0624 and 
.0832 grams if all six thrusters are operating to maintain 
attitude in their respective limit regions. The yaw and roll- 
pitch control cycle duration for this example are 108 and 72 
seconds, respectively. Assuming that 2/3 of the original 
32.35 kg of hydrazine is available for limit cycle attitude 
control, total limit cycle lifetime would be about 320 days. 
Approximate limit cycle lifetimes for various combinations of 
weighting factors, rate limits and minimum impulse bits based 
on the hydrazine budget are listed in Table 5-1. Increasing 
the value of ìà or decreasing the rate error limit will 
decrease the rate error and thus lengthen the thruster ofp 
time, but the maximum effective value of A is limited by the 


the minimum impulse bit on-period. For the thruster used 
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Table 5-1 LIMIT CYCLE LIFETIMES 


u, AT, (sec) rate(°/sec) A nit, (seem ini days) 





n- At. (Sec) rate( /sec) A SEO) 





in this example, A>100 has negligible effect since the minimum 


AV imparted is 


AV, = u, At = .0774(.025) = .001935 rad/sec 


By SC = .1172(.025) = .00293 rad/sec 


which will kick the rate error from the negative limit of 

=8.73x10" rad/sec to a positive value of .0012 rad/sec. The 
example illustrated in the figures uses an integration step 
size of .01 seconds. In order to implement the fuel saving 
choices listed in table 5-1 attitude control thrusters with 


smaller minimum impulse bits and\or on-times must be utilized. 


79 


Limit cycle lifetime based on the number of available 
pulses ranges from 30 to 90 days. This is considerably less 
than limit cycle lifetimes based on fuel consumption. In 
addition to smaller impulse bits, the number of available 
pulses must increase to lengthen the active life of ORION. 
Another option would be to operate for specified periods 
during each orbit vice continuously. Rather than replacing 
the thrusters with more complex and expensive items, this 
could increase the active lifetime of the existing design to 


acceptable levels. 


E. ENVIRONMENTAL DISTURBANCE EFFECTS 

Since the majority of the satellite lifetime will be 
spent in the operational or steady state mode, this is the 
phase when the effects of the external disturbances of gravity 
gradient and aerodynamic drag will play the largest part. 
Substituting the small angle approximation for the angular 
velocity vector and adding a disturbance moment to the Euler 


Moment Equations yields 


ó 7 a TEE 6,4 MAD (5.3999 
9. — UD 
y =u t D: 
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Il 


where a Cat Lee lp eee oc eee dts Curbanee contribution 
and u, = T./I, is the control offering. The small angle 


approximations for the gravity gradient moment equations are 


: | $ (I; 2 
To soe 9 (1,-1,) (5.34) 
0 


Enero the orbital rate w, = u/R. Figure 5-20 shovws the 
orbital rate versus altitude for any satellite orbiting the 


Earth. 
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Figure 5-20 Orbital Rate, w,, versus 
altitude 


Substituting (5.34) into (5.33) and eliminating the control 


term gives 
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, gi 
$ = -2u, | | $ (5.35) 


ys = 2 | T^ | 
0 = 30, T 0 





which are equations of motion for simple harmonic oscillators. 
The moment of inertia ratios for ORION equal +.5929 and 

-.5929, respectively. Figures 5-21 and 5-22 illustrate the 
harmonic motion about the roll and pitch axes. Plots in 5- 
21a and 5-22a are phase plane representations while graphics 
in 5-21b and 5-22b show the time response for angle and rate. 
The motion is stable and periodic with a period that 
corresponds to the satellite orbit. The orbital period is 


defined as [4] 


Rp 





(5.36) 


e 


(1) 


Figure 5-23 furnishes the orbital period versus altitude for 
an Earth orbiting satellite. In Figures 5-21 and 5-22 the 
maximum rate error induced never exceeds the initial value, 
but the angle error will exceed the desired pointing accuracy 
over the course of an orbit if a rate error exists. Gravity 
gradient torque has a time constant that is equivalent to an 


Orbital period. The time constant of this all-thruster 
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Figure 5-21a Gravitational Roll 
Effects Phase Plane 
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Figure 5-21b Gravitational Roll Effects 
Time Response 
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GRAVITY PITCH DISTURBANCE 
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Figure 5-22a Gravitational Pitch 
Effects Phase Plane 
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Figure 5-22b  Gravitational Pitch Effects 
Time Response 
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Figure 5-23 Orbital Period versus Altitude 


attitude control system is on the order of 100 seconds. If 
the control system is operating continuously to maintain the 
desired pointing accuracy, then the effects of gravity are 
insignificant compared to the normal cycling that occurs in 
the limit region. If there is zero rate error and the angle 
error is in the limit region, the gravity gradient torque will 
aid in maintaining the attitude in the limit region. 


The disturbance torques contributed by aerodynamic drag 


are 


-F, L sinécosy 
ar > 0 (5.37) 
F) L cosfcosy 
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where F, is the drag force as a function of atmospheric 
density. Table 5-2 lists atmospheric densities at some 
altitudes. [8] 


Table 5-2 ATMOSPHERIC DENSITIES 


Altitude (km) Density (kg/m?) 





At h=400 km the drag force is 2.52x10" Newtons. For small 


angles (5.37) reduces to 


pee -1.92x10° 9 N-m 


1.92x10” N-m 


Tò: = Fp D 


where L=.076m is the center of pressure offset from the center 
of mass. The effects of aerodynamic drag alone for an orbital 
period are illustrated in Figures 5-24 and 5-25 about the roll 
and yaw axes for initial conditions of angle=.01 rad and 
rate=.0001 rad/sec. The pitch angle error was assumed 
constant at .01 rad. The satellite response to drag is 
unstable, yet like gravity gradient, drag causes a slow 
response compared to the continuously operating all-thruster 
control system. Drag torques that would cause the attitude 
of the satellite to be unstable over the course of an orbital 


period do not affect motion during continuous control. 
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AERODYNAMIC DRAG ROLL EFFECTS 
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Figure 5-24a Drag Roll Effects Phase Plane 
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Figure 5-24b Drag Roll Effects Time 
Response 
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AERODYNAMIC DRAG YAW EFFECTS 
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Figure 5-25a Drag Yaw Effects Phase Plane 
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Figure 5-25b Drag Yaw Effects Time Response 
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Finally the effects of gravity gradient and aerodynamic 
drag were combined. Substituting the disturbance torques into 


(5.33) yields 


$ = -a wo $ + u, - 2.84x10° 9 (5.38) 


which is the equation of a marginally stable oscillator. The 
effect that gravity and drag have on roll is illustrated in 
Figures 5-26 and 5-27 for initial conditions of angle=.01 rad, 
zero rate and rate=.0001 rad/sec, respectively. The 
oscillation does not increase without bound and actually never 
repeats the same cycle. The phase plane plots have a run time 
of 150 minutes while the time response plots have a run time 
of 300 minutes. Again, the time constant of the gravity-drag 
torque is orders of magnitude bigger than that of the all- 
thruster control system. If the control system were only 
operated at discrete intervals during each orbit the external 
Gisturbance torques could together run the rates up to 


unacceptable levels. 
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Figure 5-26a Gravity-Drag Roll Effects with 
Zero Rate Error Phase Plane 
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Figure 5-26b Gravity-Drag Roll Effects with 
Zero Rate Error Time Response 
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GRAVITY-DRAG ROLL EFFECTS 
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Figure 5-27a Gravity-Drag Roll Effects with 
Rate Error Phase Plane 
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Figure 5-27b Gravity-Drag Roll Effects with 
Rate Error Time Response 
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VI. SUMMARY AND CONCLUSIONS 


The purpose of this thesis was to design an all-thruster 
three-axis attitude control system for ORION, the Naval 
Postgraduate School satellite bus. Pontryagin's Minimization 
Principle was employed in the derivation of fuel-time 
Minimization control laws for the acquisition and on-station 
phases of flight. The original design was modified slightly 
to accommodate a three-axis strategy. 

The original design history and specifications are 
discussed in Chapter II. The total mass of the ORION 
satellite is 250 lbm including 71.5 lbm of the monopropellant 
hydrazine for maneuvering. The attitude control system 
consists primarily of one 5 lbf thruster for orbital transfer 
and six .1 1bf thrusters for rotational control. The original 
mass placement and principal moments of inertia were left 
unchanged for simulation purposes. The four stabilizing booms 
are not utilized in this analysis. 

Additional concerns that affect a three-axis stabilized 
satellite compared to the original spin stabilized design are 
addressed in Chapter III. To meet power requirements, solar 
panels will be required vice just body-mounted solar cells. 
A three-axis ORION will be less fuel efficient than a spinning 


version. Various schemes using different combinations of 
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actuators are presented for future study. Alternate booster 
selections that could transfer ORION to higher orbits than the 
space shuttle and minimize the use of the 40 1bf thruster are 
suggested. 

The theory behind satellite motion, attitude 
determination and control hardware for a three-axis stabilized 
spacecraft, 1s developed in Chapter IV. Rigid-body spacecraft 
dynamics and the effects of the environmental disturbances of 
gravity gradient and aerodynamic drag are covered. Attitude 
control sensors and actuators are described with an emphasis 
on the higher degree of complexity required of these 
instruments vice those used on a spinning satellite. Sensors 
include horizon sensors, sun sensors, magnetometers, and rate 
gyros while the actuators are strictly thrusters. 

The derivation and simulation of the control laws that 
describe the acquisition and on-station phases of flight is 
found in Chapter V. It was determined that a dual-axis 
minimum time control law was optimal for reducing the rates 
as quickly as possible in the acquisition phase. Once the 
rates are reduced to acceptable levels, a fuel-time control 
law is implemented that utilizes a weighting factor to 
minimize fuel usage. 

Two regions of interest developed during the analysis of 
the fuel-time control law. The boundary region is a 


rectangular area centered about the desired position. The 
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deadzone is determined by the nominal pointing accuracy (+1° 
ORION design goal) in angle and the minimum thruster impulse 
bit for rate. The boundary region is distinguished by a four 
pulse limit cycle that repeats itself every 50-200 seconds for 
ORION. Assuming that 2/3 of the total hydrazine budget is 
available for limit cycle maintenance, lifetimes of 4-28 
months were obtained with all thrusters operating. The 
limiting factor for the selected thrusters was the number of 
pulses available before the Shell 405 catalyst lost 
effectiveness. The limit cycle lifetime decreased to 30-90 
days. The second region of interest was outside the boundary 
region, but close enough for the satellite to still experience 
independent or uncoupled motion about one or more axes. A 
double switching curve control law distinguished by "bang-off- 
bang" control was realized successfully. 

The environmental disturbances of gravity gradient and 
aerodynamic drag were also simulated in Chapter V. The 
gravity gradient disturbance imitates a simple harmonic 
oscillator with a period corresponding to an orbital period. 
A typical orbital period for a low earth orbit is on the order 
of 90 minutes. Aerodynamic drag applies an unstabe torque, 
but when coupled with gravity causes a random but marginally 
stable type of oscillation. Both disturbance torques are so 
small over the course of a limit control cycle, that the 


effects are negligible. 
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The three-axis stabilized design simulated in this thesis 
could be a viable attitude control option for the general 
purpose ORION bus. Upgraded thrusters must be selected that 
have smaller minimum impulse bit on-times and thrust values 
in addition to an increased catalyst lifetime if ORION is to 
operate in the constant update control mode for a feasible 
duration. Another possibility that is only mentioned in this 
work, concerns intermittant operation. ORION could function 
for specified periods every orbit when the desired target is 
in the field of view. This method would insure an extended 
lifetime, but the influence of external torques could be a 


limiting factor. 
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APPENDIX 
ATTITUDE CONTROL AND EXTERNAL DISTURBANCE PROGRAMS 


The following programs simulate a three-axis stabilized 
satellite in the acquisition and on-station phases. The 
effects of gravity gradient and aerodynamic drag are modeled 
in separate routines. 

The software package TUTSIM" was used to execute the 
control laws derived in Chapter V [11]. TUTSIM™ models 
continuous time control systems in a signal flow format. The 


following commands are defined for the potential user: 


INT - integrator block 

CON - constant block 

GAI - gain block 

MUL - multiplier block 

DIV - divide block 

ABS - absolute value block 
AND - logical and block 
OR = logical or block 

IFE - If-then statement 
REL - relay block 


The numbers preceding commands are simply block identifiers. 
The numbers following commands are block inputs. The constant 
block does not have an input. The integrator and gain blocks 
require at least one input. The multiply, divide and logical 
operator blocks need at least two inputs. The if-then 
statement has three inputs. The first input is the argument. 


The second and third inputs are the commands that support true 
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or false arguments, respectively. The relay block has four 
inputs. The last input is the argument. The first three 
inputs are commands that support greater than, equal to or 
less than results, respectively. Further information on the 
capabilities of TUTSIM™ can be found in the operator manual. 

The first program simulates on-station control about one 
axis. Input parameters include initial conditions, principal 
moments of inertia, control torques, limit region boundary 
values and weighting factors for two switching curves. The 
command logic applies control according to the satellite 
attitude in the phase plane in relation to the switching 
curves and boundary limits. The control law for this program 
is derived in the first section of Chapter V. 

The next three programs simulate the acquisition phase 
for single and double axis control. Input parameters include 
initial conditions for rate errors, principal moments of 
inertia and torque values. These programs apply control as 
a function of the rates for coupled equations of motion. 
Switching occurs about the ordinate and abscissa in a defined 
outer region. Switching occurs on minimum time curves in the 
inner region. The control law for the acquisition routines 
can be found in the second section of Chapter V. 

The last two programs simulate the effects of gravity 
gradient and aerodynamic drag distubances, respectively. 


Inputs parmeters are provided for initial conditions, 
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principal moments of inertia and orbital altitude. The last 


section of Chapter V illustrates the effects of external 


disturbances. 


omony nn & Un 


33 


34 
35 
36 
37 
38 
39 
40 


41 
42 
43 
44 


INT 
INT 
CON 
CON 
DIV 
CON 
CON 
CON 
ABS 
ABS 
CON 
CON 
CON 
CON 
MUL 
MUL 
SUM 
SUM 


AND 
MUL 
MUL 
MUL 
CON 
IFE 
IFE 


IFE 
IFE 
IFE 
IFE 


n 
WO 0 “J 009 00 


29 
25 
26 
27 


34 
28 


29 
37 
36 
BD 


mn 


31 


31 


40 
41 


43 
44 
44 
39 


:X1,ANGLE 
;X2, ANGLE 
SIX Iy, or 
SIX ey OL 
;Ux,Uy,or 
a 

sb 

“4 

; ABSOLUTE 
: ABSOLUTE 


RATE 
Iz 
Tz 
Uz 


ANGLE 
RATE 


ANGLE LIMIT 
RATE LIMIT 


: LAMDA1 
; LAMDA2 
4% LAMDAl 
; G*LAMDA2 


: 1+4 LAMDA1 
: 144 LAMDA 2 


oa 


‚.5(1+4LAMDA1)/U 
: . S(1*4 LAMDA2) /U 


: 9x2 [X2] 


:A1*X2[X2] 
:A2*X2[X2] 


:X14,5X2[ X21 
: X14-A1*X2(X2] 
:X1+A2*X2(X2] 


:X1*Xx2 
;X1«ANGLE 


LEMEE 


;X2«RATE LIMIT 
;-Ux, -Uy or -Uz 
; U*SIGN (X1) 


;U*SIGN(X2) 


:X1«ANG LIM AND X2«RATE LIM 


: ABOVE OR 
: ABOVE OR 
: ABOVE OR 
;U-0 

SATTITUDE 
"ATTITIUDE 


;X1<ANGLE 
: RELATION 
: RELATION 
: RELATION 


1 On-Station Control 
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BELOW MIN TIME CURVE 
BELOW 2nd CURVE 
BELOW 3rd CURVE 


IN LIMIT REGION? 
IN 1st/3rd QUADRANT? 


LIMIT? 

TO 3rd CURVE? 

TO 2nd CURVE? 

TO MIN TIME CURVE? 


1 INT 3 ;Wx 
2 INT 4 27 Wz 
SMU 17 2 8 ;(Iy-Iz)/Ix*Wy*Wz 
4 MUL 18 I 8  ;(Ix-Iy)/Iz*Wy*Wx 
5 ABS 1 : ABSOLUTE ROLL RATE 
6 ABS 2 : ABSOLUTE YAW RATE 
7 CON i 
8 CON ¡Wy=.01 
9 CON "TX 
10 CON SI 
11 CON “12 
12 CON Te 
13 DIV 12 11 AE 
14 GAI 13 ;-Uz 
15 SUM 10 -11 ‚Iy-Iz 
16 SUM 9 -10 ‚Ix-Iy 
17 DIV 15 9 ;(Iy-Iz)/Ix 
18 DIV 16 LI ,(Ix-Iy)/Iz 
19 MUL 7 2 6 ;.5Wz[Wz] 
20 MUL -7 i 5 ;-.5Wx[Wx] 
21 SUM 1 20 19  ;Wx-.5Vx[Vx]*.5Wz[Wz] 
22 MUL 1 I :Wx^2 
23 MUL 2 2 :Wz^2 
24 MUL 29 13 13 ;4*Uz*2 
25 SUM 22 23 -24 ;Wx*24+Wz*2>4*Uz*2 
26 REL 14 13 13 ;U=-U*SIGN(Wx- . 5Wx[Wx]+.5Wz[Wz]}) 
21 
2/<TEE 25 28 26 ¡INNER CONTROL 
28 IFE 2 14 13 ;OUTER CONTROL 
29 CON “4 


2 Acquistion Control for a Single 
Actuator Set, Uz 


29 


1 INT 3 27 Wx 
2 INT 4 ;Wz 
3 MUL 17 2 8. ;(Iy-Iz)/Ix*Wy*Wz 
4 MUL 18 1 8  ;(Ix-ly)/1z*Wy*Wx 
5 ABS 1 : ABSOLUTE ROLL RATE 
6 ABS 2 : ABSOLUTE YAW RATE 
7 CON P ur 
8 CON SMy-.01 
9 CON a 
10 CON 35" 
11 CON I2 
12 CON EX 
13 DIV 12 9 ;Ux 
145 GAI 13 : -Ux 
15 SUM 10 11 ‚Iy-Iz 
15 SUM 9 -10 ‚Ix-Iy 
17 DIV 15 9 ; (Iy-Iz)/Ix 
18 DIV 16 11 ¡(LX TY) / TZ 
19 MUL 7 2 6 2, SWZ (We) 
20 MUL 7 1 5 ;.5Wx[Wx] 
21 SUM -2 19 20 ;-Wz+.5Wx[Wx]+.5Wz[Wz] 
22 MUL 1 1 :Wx^2 
23 MUL 2 2 “Wz 2 
24 MUL 29 13 13 :4*Ux^2 
25 SUM 22 23 -24  ;Wx^24Wz^2»4*Ux^2 
26 REL 14 13 13 ;U--U*SIGN(-Wz*.5Wx[Wx])4.5Wz[Vz 
21 
2I IFE 25 28 26 ¡INNER CONTROL 
28 IFE I 14 13 ;OUTER CONTROL 
29 CON :4 


3 Acquisition Control for a Single 
Actuator Set, Ux 
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26 REE 


Z TREL 


28 MUL 
29 MUL 
30 CON 
JSISMUE 
32 ESUM 
OR EL 


34 REL 


SS IFE 
36 IRE 


4 Acquisition 


Actuator 


14 
28 
17 


15 


32 
32 


35 -Wx 
36 -Wz 
2 8 ;(Iy-Iz)/Ix*Wy*Wz 
E 8 ;(Ix-Iy)/Iz*Wy*WUx 
;ABSOLUTE ROLL RATE 
¡ABSOLUTE YAW RATE 
ND 
;Wy=.01 
;Ix 
‚Iy 
-Iz 
S 
SEZ 
9 ;Ux 
; -Ux 
11 "Uz 
‚-Uz 
-11 ‚Iy-lz 
-10 V IxeIy 

9 ; (Iy-Iz)/Ix 
LI very Liz 

1 5 Wx (Wx } 

2 6 ;.5Wz[Wz] 

2 -22 ¡Wx+Wz-.5Wx[Wx)]+.5Wz[Wz] 
-2 22 ;Wx-Wz+.5Wx[Wx]+.5Wz[Wz] 
16 16 ;-Uz*SIGN(25) 
14 14 ;-Ux*SIGN(24) 

1 ‚Wx”2 

2 2200 

7-4 

14 30 ;-4*U^2 
29 31 ;Wx^24Wz^2»54*U^2 
16 16 ;Uz--Uz*SIGN(Wx) 
14 14  ;Ux--Ux*SIGN(Wz) 
33 26" UZ 
34 27 Ux 

Control for a 
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Double 


INT 
INT 
INI 
INE 
CON 
CON 
CON 
CON 
CON 
10 CON 
11 SUM 
12 CON 
13 CON 
14 SUM 
15 SUM 
16 DIV 
17 DIV 
16 MUL 
19 DIV 
20 MUL 


DW 0 -JOwvUn £^ UL F9 


21 MUL 
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l4 


11 


12 


13 


10 
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-5 


11 


19 


19 


11 


16 


17 


POT 

,D(PHI)/DT 

; THETA 

; D(C THETA) /DT 

ra 

ly 

"TI 

;GRAVITATIONAL PARAMETER-3.986e5 
;HEIGHT 

‚EARTH RADIUS=6378 KM 
;‚R=Re+H 

2-2 

i3 

‚iv Iz 

CISCIX 

,(Iy-Iz)/Ix 

; (Iz - IX)/Iy 

R3 

;ORBITAL RATE,Wo^2-MU/R^3 
“2 Nom? (Iy-Iz)/Ix*PHI 


‚3 Wo^2 (Iz-Ix)/Iy*THETA 


5 Gravity Gradient Disturbance Effects 
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INT 
INT 
INT 
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20 
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22 
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; PHI 
; D(PHI)/DT 
: PSI 
;D(PSI)/DT 
TX 
25/ 
“EZ 
:GRAVITATIONAL PARAMETER=3.986e5 
"EARTH RADIUS=6378 KM 
: HEIGHT 
:R=Re+H 
‚Iy-Iz 
; (Iy-Iz)/Ix 
;DENSITY,RHO 
:ORION CROSS-SECTION,A-.429m^2 
ied] 
;DRAG COEFFICIENT-2 
;ORBITAL VELOCITY,V^2 
30 ;DRAG FORCE 


:L-.076 m 
DIOS RS 
;ORBITAL RATE,Wo^2 
:DRAG YAW TORQUE 
: THETA 
:DRAG ROLL TORQUE 
1 ;Wo^2 (Iy-Iz)/Ix*PHI 
:FD*L/IZ 
;FD*L*THETA/Ix 
;le6 
U 2TIN M 2/572 


Aerodynamic Drag Disturbance Effects 
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